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Abstract 

Discrete nonlinear Schroginger equation (DNLS) of the form, z^p- = C n +\ 

+ C n -\ - Xn[\C n +i\ 2 + |C„_i| 2 — 2|C n | 2 ]C n is used to study the formation of 

stationary localized states in one dimensional system due to a single as well 

as a dimeric nonlinear impurity. The fully nonlinear chain is also considered. 

The stability of the states and its connection with the nonlinear strength is 

presented. Results are compared with those obtained from other DNLS. It is 

found that the DNLS used in this paper has more impact in the formation of 

stationary localized states. 

PACS numbers : 71.55.-i, 72.10.Fk 
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I. INTRODUCTION 



It is well know that the transport properties of a system is directly related with the 
formation of localized states in the system. Localized states appear due to the presence of 
impurity or disorder (which breaks the translational symmetry) in the system ffl. There 
has been studies on the formation of localized states due to linear impurities in various 
systems ||T| . On the other hand, only a few look into the formation of localized states due 
to nonlinear impurities. The discrete nonlinear Schrodinger equation, used to study the 
formation of stationary localized (SL) states P-|T3|] is given by 

.dC„ 



dt 



£nC n + V(C n+ i + C n ^\) — Xn\C n \ a C n , (1) 



here C n is the probability amplitude of the particle (exciton) to be at the site n, e n and 
Xn are the static site energy and the nonlinear strength at site n respectively and V is the 
nearest neighbor hopping element. The nonlinear term \C^CC n arises due to the interaction 
of the exciton with the lattice vibration |14,15||. The above eq. (1), has been used to study 



the formation of SL states in one dimensional chain as well as in Cayley tree with single 
and dimeric nonlinear impurities p|-|6|,p|-p!l|1 . For the case of a perfectly nonlinear chain 
where \ n = x, it has been shown that SL states are possible even though the translational 
symmetry of the system is preserved |7j,P|,|n|-|TB|] . These results were remarkably different 
when compared with that of the corresponding systems with linear impurities. 

The equation (1) has been derived with the assumption that the lattice oscillators in the 
system are local and oscillate independently. A natural question to ask is, what happens to 
the formation of SL states when the oscillators in the lattice are occupied with their nearest 
neighbors. In this case the discrete nonlinear Schrodinger equation takes the form 

4C n 



dt 



£nC n + + Cn— l) — Xn(|C n +i| + |C n _i| — 2\C n \ 2 )C n (2) 



where C n , e n , V and Xn carries the same meaning as in eq. (1). We notice that the eq. 
(2) has more nonlinear terms compared to the eq. (1). To the best of our knowledge, this 



2 



equation has not been used to study the formation of SL states even though eq. (2) is more 
important in condensed matter physics. Our intention is to look for the formation of SL 
states in one dimensional system due to the presence of nonlinear impurities (as described 
by eq. (2)) and further to compare the results with those obtained from eq. (1) and to see 
which one has more impact in the formation of SL states. 

The organization of the paper is as follows. In sec. II we discuss the effect on the 
formation of SL states due to a single impurity (i.e. Xn = x(^n,o))- 111 sec - HI we consider 
the case of dimeric impurity (i.e. Xn = x(^n,o + <W)) anc ^ * n sec ^ we consider the perfectly 
nonlinear chain. In sec V we discuss about the stability of the SL states. Finally in sec. VI 
we summarize our findings. 

II. SINGLE NONLINEAR IMPURITY 

Consider the system of a one dimensional chain with a nonlinear impurity at the central 
site. The time evolution of an exciton in the system is governed by eq.(2) with Xn = X^n,o- 
The Hamiltonian which can produce the equation of motion for the exciton in the system is 
given by 

H = £(<?*C n+ i + C n C* n+1 ) - I [Idl 2 + IC^I 2 - 2|Co| 2 ]|C | 2 . (3) 

As J2n \C n \ 2 is a constant of motion, we suitably renormalized so that J2n |Cn| 2 =L We call 
it normalization constant. Therefore, \C n \ 2 can be treated as the probability for the exciton 
to be at site n. Since we are interested in finding the stationary localized states, we consider 
the ansatz 

C n = <j> n exp(-iEt); <p n = <M H (4) 

where < r) < 1 and i] can be asymptotically defined as r] = \ E \~^ E ' 2 -i ^ jt; [ s the energy of 
the localized state which appears outside the host band. Since in a one dimensional system, 
states appearing out side the host band are exponentially localized, the ansatz (in eq. (4)) is 
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justified as can also be readily derived from the Greens function analysis PJlO|. Substituting 
the ansatz in the normalization condition we get 



2_ W 

" 1 + rf ' 



(5) 



Direct substitution for O , <p n and hence C n in terms of r\ in eq. (3), yields an effective 
Hamiltonian, 

H *ff ~ TT^f + (T+W ' ( } 

The fixed point solutions of the reduced dynamical system described by H e ff will give the 
values of rj (which correspond to the localized state solutions) |§. Note that the effective 
Hamiltonian is a function of only one dynamical variable, namely, rj as x is constant. Thus 
fixed point solutions are readily obtained from the condition dH e ff/dr] = 0, i.e., 

4 = .KWHio + V) = 

X 0- + T) 

Thus the different values of r\ G [0,1] satisfying the eq. (7) will give the possible SL states 
for a given value of \. It is clear from the expression for f(rj) that f(t]) — > as r/ — > 
and 77 — 1 . Therefore it is expected that /(^) will have at least one maximum, which is 
indeed the case as can be seen from Fig. (1) where f(rj) is plotted as a function of 77. Notice 
that there will be no graphical solution if ^ > f(r) max ), one solution if ^ = f(r) max ) and 
two solutions if - < f(r] max ). Thus there is a critical value of x> say, Xcr below which no 
localized states are possible and is given by 

x " = JOTAJ = 1 ' 2696 - (8) 

Thus for \= 1.2696, we get one SL state and two for \ > 1.2696. For a system described 
by eq. (1) (with a— 2), it has been shown in ref. that the corresponding critical value 
for x is 2. Also the maximum number of SL states possible was 1. Thus we see that the 
nonlinearity arising in eq. (2) reduces the critical strength and produces more number of 
SL states. Hence, the eq. (2) is indeed more effective in the formation of SL states. 



III. DIMERIC NONLINEAR IMPURITY 



We consider the case where the one dimensional lattice has two nonlinear impurities at 
site and 1 respectively, i.e., x(^n,o + <5n,i)- As in the case of single impurity, it is easily 
verified that the Hamiltonian for the system is given by eq. (2) with x as defined above. 
For stationarity condition we assume that C n = <f) x exp(—iEt). Furthermore, for localized 
states we assume the following form for <f> n . 

n = [sgn{E)rj\ n ~ l <j)\] n > 1 

and 

0_ w = [sgn(E)r]}^(f) ; n<0 (9) 

with 1] as defined earlier. The ansatz is justified as those states which appear outside the 
host band are exponentially localized (which can be derived exactly from the Greens function 
analysis ||). Three different possibilities arise, (i) 4>i — 0o (symmetric case), (ii) 4>i — — <Po 
(antisymmetric case) and (iii) <pi ^ <po (asymmetric case). It is possible to encompass all 
the different cases by introducing a variable /3 = 4^. The value of (3 is confined between 1 
and -1 if |0o| < 4>i- Else we inverse the definition of (3. (3 = 1, —1, and ^ 1 correspond to the 
symmetric, antisymmetric and the asymmetric state respectively. Substituting the ansatz 
as well as the definition of (3 in the normalization condition, J2°?oo \Cn\ 2 = 1 we get 

W = ^J do) 

and the reduced Hamiltonian 
If (3 = ±1 we get, 

Hf ff = T 2 V + 2sgn(E) + 3 X± r/ (1 ~ ^ (12) 
Here '+' sign corresponds to the symmetric case and '-' sign corresponds to the antisymmetric 



case. The number of fixed point solutions of the reduced dynamical system described by H e ff 
gives the the possible number of SL states. The fixed point solutions satisfy the equation, 



— = 377(1 =F77)(1 ±V 2 )- (13) 
X± 

From eq.(13) it is clear that there exists two critical values of \ namely, 0.7149 and 1.6525. 
There is no SL state for x < 0.7149, one symmetric SL state at x — 0.7149, two symmetric 
SL states for 0.7149 < x < 1-6525, two symmetric and one antisymmetric SL state at % = 
1.6525 and two symmetric and two antisymmetric SL states for x > 1.6525. 

Now let us consider asymmetric case where (3^1. The effective Hamiltonian is function 
of two dynamical variables namely, (3 and rj. Therefore the fixed point solutions will obey 
the equations given by 

OB-ff 0and^£ = 0. (14) 



dr\ " dp 

After a little algebra we obtain the desired equation, 

1 (3(9 - 7(3 2 -(3 4 - (3 G ) 2 



/(/?)■ (15) 



X 2(1 -/P)(3 -/*»)< 
The function f((3) monotonically increases with (3 and it goes to infinity as (3 goes to 1. 
From this we can immediately see that there always exists one SL state no matter how small 
X may be. 

Combining all the possible states we find that there is one SL state for x < 0.7149, two 
at x = 0.7149, three for 0.7149 < x < 1-652, four at x = 1-6525 and five for x > 1-6525. 
Hence the maximum number of SL states is five. We further note that the critical values for 
nonlinear strength is lower and the number of SL states are more compared to the results 
obtained by eq. (1) (with a— 2) |L0|| . Thus it is again confirmed that eq. (2) is more effective 
in the formation of SL states compared to eq. (1). 



IV. FULLY NONLINEAR CHAIN 

We now consider perfectly nonlinear chain i.e., Xn = X- The Hamiltonian for this 
system is given by eq.(2) with Xn = X- Using the stationarity condition, we can obtain 
the Hamiltonian in terms of (b n . In this case it is not possible to find the exact ansatz for 
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the localized states, but there are a few rational choices. For example, a single site peaked 
as well as inter-site peaked and dipped solutions are possible. We will consider these cases 
subsequently. Let us first consider the on-site peaked solution. Without any loss of generality 
we can assume that the exciton profile is peaked at the central site. Therefore, using the 
ansatz n = 0O 7 ?'™' an d the normalization condition we get the effective Hamiltonian, 



From the fixed points equation, dH e ff/dr] = 0, we obtain 

1 3r/(l-r/ 2 ) 



(17) 



X (1 + r/ 2 ) 2 

After analyzing this equation we find that there is a critical value of x — 1-333 below 
which there is no SL state and above it there are two states. At the critical value of x there 
is one state. 

For the inter-site peaked and dipped solutions we use the ansatz of the dimeric impurity 
nonlinear impurity, carrying out the calculation involved we obtain the effective Hamiltonian 
of the reduced dynamical system to be 

where (3 is defined earlier. We first consider the case (3 = ±1. Substituting f3 — ±1 into the 
Hamiltonian and from the fixed points equations we obtain 

J_ = r / (l-r / 2 ) 2 (2 + r/ 2 ) 
X ± 2(sgn(E)Tv)(l + V 2 ) 2 ' 

Here '+' sign corresponds to the symmetric case and the '-' sing to the antisymmetric case. 

From eq. (19) it is clear that there will be two critical values of x namely, xtr = 2.4653 and 

X„ = 5.9178. There is no SL state for x < xtn one SL state for x — xtn ^ wo ^ states for 

xtr < X < Xcn three SL states at x — Xcr an< i f° ur SL states for x > Xcr- 

On the other hand for (3 ^ 1 we find that (3 G [0, 1] and r\ e [0, 1] satisfy the following 

equations. 



- 4/^(1 + /3 2 )(1 + rff + 2sgn{E){\ + /? 2 ) 2 (1 + rff + 2 X?7 [-3 + (3 2 - f3 2 V 4 } 

- X [2r] 2 - 4/3V - 2/3 4 + r? 4 - 2/3V] = 
2(1 + r ? 2 )(l + P 2 ) 2 - 4/3 2 (l + r ? 2 )(l + /? 2 ) - x [(l + /3 2 )(2/3 - 2p V tfr? - 4/5 3 ] 

+4x/9[/9 2 - /5 2 r/ 4 + rf + 0V - 1 - (3 4 ] = (20) 

As it is not possible to decouple the equations, we have obtained numerically the possible 
values of (5 and rj for various values of x- It is found that there is always exists one SL state 
for any nonzero value of x- 

Now combining all the possibilities, we obtain the following result for the fully nonlinear 
chain. There will be only one SL states for x < 2.4653, two for x — 2.4653, three for 
2.4653 < x < 5.9178, four for x = 5.9178 and five for x > 5.9178. Hence the maximum 
number of SL states is five. We further note that SL state appears even if the system is 
perfect (the translational symmetry is preserved). Therefore, we may call these states to be 
self-localized states. 



V. STABILITY 

The stability of the SL states can be understood from a simple graphical analysis. For 
this purpose, consider the case of single impurity with x = 1.3 (for which two SL states 
appear). The fixed point equation for the single impurity case with x — 1-3 is given by 

= i + = . (21) 

(1 + i] z ) 

The flow diagram of the dynamical system described by the H e ff given in eq. (6) is con- 
structed in the following manner. We treat G(rj) as the velocity and rj as the coordinate of 
the dynamical system. G{rj) is plotted as a function of rj in fig. (2). 'A' and 'B' are the 
fixed points corresponding to the SL states. If G(rf) > 0, the flow of the dynamical variable 
is in right direction else it is in left. The direction of flow is shown by arrows in different 
regions. It is clear from the flow diagram 'A' is a stable fixed point whereas 'B' is unstable. 
Therefore, for the case of single impurity, one state is stable and the other one is unstable. 
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The energy of the SL states as a function of \ is plotted in fig. (3). Once again we 
confine to the single impurity case. It is clear that energy of one state increases and that of 
the other decreases. (Note that the points 'A' and 'B' of fig. (3) gets mapped to points 'A" 
and 'f?" respectively.) Thus we conclude that the states in the upper branch of the energy 
diagram are stable SL states and those of the lower branch are unstable SL states. In other 
words, if the energy of SL state increases with the increase of nonlinear strength, the state 
is stable otherwise, unstable. 

VI. CONCLUSION 

DNLS given by eq. (2) is used to study the formation of stationary localized states in 
a one dimensional system with single and a dimeric nonlinear impurity. It is found that 
the number of SL states are more than the number of impurities in the system. Maximum 
number of SL states due to single nonlinear impurity is two and that due to dimeric nonlinear 
impurity is five. It is further found that SL states may appear even in a perfectly nonlinear 
system. Thus one may call these SL states as self- localized states. It is also interesting to 
note that eq. (2) is more effective in the formation of SL states compared to eq. (1). The 
stability of the SL states is discussed and the connection of the stability of a state with its 
energy variation as a function of nonlinear strength is presented. For a clearer understanding 
on the effect of nonlinear impurities on the formation of SL states, one needs to consider the 
presence of a finite nonlinear clusters in a linear host lattice. Investigation in this aspect is 
in progress and will be reported elsewhere. 
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FIGURES 



FIG. 1. f(rj) is plotted as a function of r\. It clearly shows a maximum. 

FIG. 2. G(rj) is plotted as a function of ry. ^4 and are fixed points. 
FIG. 3. Energies of the SL states are plotted as a function of %. There is no state below Xcr 

and two states above it. Energy of one state increases and that of the other decreases. Points A' 

and B' corresponds to points A and B respectively of Fig. 2. 
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